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Abstract 

We construct a class of topological excitations of a mean field in a two- 
dimensional spin system represented by a quantum Heisenberg model with 
high powers of exchange interaction. The quantum model is associated 
with a classical one (the continuous classical analogue) that is based on 
a Landau-Lifshitz like equation, and describes large-scale fluctuations of 
the mean field. On the other hand, the classical model is a Hamiltonian 
system on a coadjoint orbit of the unitary group SU(2s + l) in the case 
of spin s. We have found a class of mean field configurations that can be 
interpreted as topological excitations, because they have fixed topological 
charges. Such excitations change their shapes and grow preserving an 
energy. 

1 Introduction 

According to Mermin and Wagner [IJ there is no ferromagnetic or antiferromag- 
netic order in the one- and two-dimensional isotropic Heisenberg models with 
interactions of finite range at nonzero temperature. This statement is proven 
due to Bogoliubov's inequality in the general case. Here we construct excitations 
that cause a destruction of a long-range nematic or mixed ferromagnetic-nematic 
order. This is an extension of the results of Belavin and Polyakov [5] • 

We model a planar magnet by a square atomic lattice with the same spin s 
at each site. We describe this two-dimensional spin system by a generalized 
Heisenberg Hamiltonian, taking into account high powers of the exchange inter- 
action (S n , S n -\s), where S n is a vector of spin operators at site n. By a mean 
field approximation we obtain a classical long-range equation from the quantum 
Heisenberg one. 

An equation for a mean field (the field of magnetization and multipole mo- 
ments) is a Hamiltonian equation on a coadjoint orbit of Lie group. At the 
same time, this is a generalization of the well-known Landau-Lifshitz equation 
for a magnetization field. In this context we obtain effective Hamiltonians for 
the magnetic system in question. Using Kahlerian structure of coadjoint orbits, 
we construct effective Hamiltonians such that their minimums are proportional 
to topological charges of excitations. In addition, we produce these mean field 
configurations that give minimums to the Hamiltonians. 
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2 Quantum and classical models 



As mentioned above, we represent the spin system by a planar atomic lattice 
with the same spin s at all sites. We assign three spin operators (S„, S^, S^) = 
S n to each atom n; the operators obey the standard commutation relations 

where a, b, c run over the values {1, 2, 3}, and 5 nm is the Kronecker symbol. 
2.1 Generalized Heisenberg Hamiltonians 

We are interested in so called high spins s > 1. In this case, we can describe the 
system by the following bilinear-biquadratic Hamiltonian 

n 2 = - ^(j(S n , S n+5 ) + K(S n , s n+5 ) 2 ) . 

Here <5 runs over the nearest-neighbour sites, n runs over all sites of the lattice, 
constants J and K denote exchange integrals. In the case of spin s^3/2, the 
above Hamiltonian can include also the bicubic exchange, namely 

ii 3 = - y^^JjSn, S n+ s) + K(S n , S n+ s) 2 + L(S n , S n+ s) 3 ^j , 

n.S 

where L denotes the corresponding exchange integral. 

One can easily write a generalized Heisenberg Hamiltonian for the system of 
an arbitrary spin s or greater than s. This Hamiltonian contains all powers of the 
exchange interaction up to 2s. It can be reduced to a bilinear form if one takes 
the 2s + 1-dimensional space of irreducible representation of the group SU(2). 
The spin operators {5°} over this space generate a complete associative matrix 
algebra, which has a sufficient number of operators to reduce the corresponding 
Hamiltonian to a bilinear form. 

For example, in the case of spin s = 1 the appropriate space of representation 
is 3-dimensional, and we choose a canonical basis in the form: {|+1), |— 1), |0)}. 
The spin operators {S^} generate the algebra Mat3 X 3- In order to form a basis 
in the algebra we take the tensor operators of weight 2 

Q°a = Sn&n + ^n^rn a 7^ &j , , 

Qf = (S 1 n f-(Slf, Ql° = V3((S*f - I) 

in addition to the spin operators. The introduced operators are called quadrupole 
operators. 

In the case of spin s = 3/2 the appropriate space of representation is 4- 
dimcnsional, and {| + |), I + \ — \), \ — §)} is a canonical basis. We complete the 
associative matrix algebra Mat4 X 4 of {S®} by the tensor operators of weights 2 
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and 3, defining them by the following formulas: 

<£* = $((#) a -(#) a ), ^° = f ($D 2 -i), (2) 

ff = (QfS* + S 3 n Qf), a, b € {1, 2}, a ± b, 

rpab 1 / / Qa\2 ob , cb ( qci\2 , na ob qci / ob \3 \ 

n — x /ft\\°n) D n~T~ D n\ D ?i) ~r °n°n°n \ D n) J> 

, \ (3) 

f " a = VTo (^M + S 3 n Qf + V3(Q 2 n °S a n + S a n Ql )) , 

f™ = ±(41S 3 n -20(S 3 n ) 3 ). 

We call the tensor operators of weight 3 sextupole operators. In what follows we 
denote all tensor operators over the chosen space of representation by {P%}. 

In terms of representation operators a generalized Heisenberg Hamiltonian 
gets a bilinear form. For the Hamiltonians considered above we have: 

n.S b 7i,5 a 

H 3 = -(J-l K +^L)J2Y,^n + s - 7 i(4K-L)N- 

n.S b 

~ I (K-2L) E ]T Q a n Q« +5 ±L Y E f % f n +S , 

n,5 a n,S /3 

where N is the overall number of sites in the lattice. Obviously, the obtained 
bilinear Hamiltonians are SU(2)-invariant, because they are constructed from 
representation operators of the group SU(2). 



2.2 Mean field approximation 

Here a mean field is a field of expectation values for the operators cal- 
culated after spontaneous breaking of symmetry. The breaking of symmetry is 
performed by switching on an external magnetic field, which specifies an order 
in the system; then the external field vanishes. Such kind of averages is also 
called quasiaverages [3]. 

We denote a mean field averaging by (•), and components of the mean field 
by {fJ*a(x n) = {P%}}. A mean field approximation of the bilinear-biquadratic 
Hamiltonian has the form 

3 8 

^MF - -{J-\K)Z Y, Y ~ \KZ Y, Y P r>a(x n ) - %KzN, 

n a—1 n a— A 

where 2 is a number of the nearest-neighbour sites. 

Evidently, a mean field Hamiltonian remains SU(2)-invariant. Then by an 
action of the group SU(2) it can be reduced to a diagonal form. Of course, 
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this reduction is possible only in the case of thermodynamical equilibrium and 
an infinite lattice, when the mean field becomes constant and the dependance 
on site n can be omitted. Moreover, almost all components of the mean field 
vanish, except the components corresponding to diagonal operators of {-P,"}- 
For the bilinear-biquadratic Hamiltonian a diagonal form is the following: 



■H 2 MF = -zN ([J-\K)S Z ^ + \KQ 20 ^ + fif) . 



The remaining components are suitable to serve as order parameters. The com- 
ponent /X3 describes a normalized magnetization (a ratio of the z-projection of 
magnetic moment to a saturation magnetization). The components /is and ^15 
are normalized projections of quadrupole and sextupole moments, respectively. 

Possible values of order parameters can be obtained from the self- consistent 
equations 

TrP a e-T^ 
Ha = (P ° MF = 1 

Tre"^ 

for all diagonal operators P a . Here we use the density matrix with the one-site 
Hamiltonian Hmf = T~Lmf /N . Note, that for the standard Heiscnbcrg Hamilto- 
nian, when only the spin operators are considered, a self-consistent equation 
turns into the well-known Weiss equation. 

In the case of bilinear-biquadratic Hamiltonian, an analysis of self-consistent 
equations gives the following. We adduce probable values of order parameters 
in the limit T^O as J, K > 0: 1) |/i 3 | = l, Ms = ; 2) |/i 3 | = |, Ms = ^^; 
3) (J.3 = 0, |/z g | = -^L; 4) ^3 = 0, |/z g | = -4=. Two of the solutions correspond 
to states with the ferromagnetic order, because they have a nonzero magnetiza- 
tion The other two solutions correspond to states with the quadrupole order 
(so called nematic states), which are states with zero magnetization. Solutions 
2) and 4) are unstable and called partly ordered. If K becomes negative, only 
solution 1) remains. These results accord with the results of [HE] and with the 
phase diagram of ordered states in a one-dimensional spin system from |6j . 

In what follows we consider an SU(3)-invariant system with the bilinear-bi- 
quadratic Hamiltonian. The SU(3)-invariance is reached by assigning J = K; 
this is the boundary line between the ferromagnetic and the nematic regions 
(see [S]). It means the system can appear in the both states. In the case of 
Hamiltonian "H 3 , the maximal SU(4)-invariance is reached as J = — 37 -£T = — %L, 
that is located within the ferromagnetic region. 

Now we apply the mean field averaging to the quantum Heisenberg equation 

ih^ = [Pz,n}. (4) 

The averaging is performed with the assumption of zero correlations between 
fluctuations of {P°} at distinct sites: (P^P^) = (P£)(P£). Then we take a 
large-scale limit and obtain the Landau-Lifshitz like equation 

— 2J7 Cabc^bijJ-c.xx + fJ-c.yy), (5) 
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where C a bc are structure constants of the Lie algebra of with the commu- 
tation relations [P°,P^] = iC a b c PnS nm . Equation Q is an equation of motion 
for the mean field {fi a (x)} over the plain {x = (x,y) | x, y£ K} that replaces 
the lattice. 

In the case of standard Heisenberg Hamiltonian ([S]) coincides with the Landau- 
Lifshitz equation for an isotropic magnet. Therefore, in the general case we call 
(O a generalized Landau- Lifshitz equation for the vector field {/i a }- The vector 
field has 8 components if one exploits the bilinear-biquadratic Hamiltonian, and 
15 components for the Hamiltonian with the bicubic exchange. 



2.3 Effective Hamiltonians on coadjoint orbits 

The generalized Landau-Lifshitz equation ((SJ can be interpreted as a Hamilto- 
nian equation on a coadjoint orbit of Lie group. In the case of spin s = 1 we deal 
with the group SU(3), in the case of an arbitrary spin s the group is SU(2s+l). 
Note, that the matrices {P a } serve as a basis in the corresponding Lie algebra 
su(2s+l), and components of the mean field {fJ, a } serve as coordinates in the 
dual space to su(2s+l). 

We start with brief description of the groups SU(3) and SU(4). For more 
material see, in particular, [5]. 

The group SU(3) has two types of orbits: the generic xj(i)^u(i) °^ dimen- 
sion 6, and the degenerate su^xufi) °^ dimension 4. Each orbit of SU(3) is 
defined by two numbers to and q, which are values of the coordinates 113 and 
fis at an initial point (a point in the positive Weyl chamber). Simultaneously, 
these numbers are limiting values of the mean field components /13 and /i 8 at 
zero temperature. For a degenerate orbit one has to assign to = 0, or to = y/3 q. 
Evidently the degenerate orbits with to = are domains of mean field configura- 
tions that realize nematic states. Ferromagnetic states are realized on all other 
orbits of SU(3). 

The group SU(4) has four types of orbits: the generic u(i)x S tj(i)xU(i) °^ 

dimension 12, the degenerate su(2)xu(ijxu(i) °^ dimension 10, the degenerate 

Sfu^xu^)) °^ dimension 8, and the maximal degenerate su^xufi) °^ dimen- 
sion 6. Each orbit is defined by numbers to, q, p, which are limiting values of the 
mean field components /13, /ig, /115. Almost all orbits are domains of ferromag- 
netic mean field configurations. Nematic states are realized on the degenerate 
orbits of dimension 8 as to = p = and q is arbitrary. So it is probable to reveal 
a nematic state even in the ferromagnetic region of the phase diagram [6 a . 

As shown above, limiting values of diagonal components of a mean field serve 
as order parameters. Simultaneously, they define a coadjoint orbit where the 
corresponding mean filed configuration lives. 

Each orbit possesses a Hamiltonian system with an equation of motion and 
a group-invariant Hamiltonian. For a degenerate orbit of SU(3) the equation is 

= 3{mi+q 2 ) CabcfJ-b(Pc,xx + f-c,yy)i (6) 
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where A denotes a dimensional constant. The values m, q of order parameters 
(a magnetization and a projection of quadrupole moment) define an orbit via 
the following equations, which we call constraints: 

SabfJ-a^b = m 2 + q 2 

7 I / m 2 +q 2 
UabcfJ-b^c — ±y 5 Ma- 

Here d aoc = Tr(P a P D P c + PbP a P c ) is a symmetric tensor. The corresponding 
SU(3)-invariant Hamiltonian is the following 

«rf - / E((Ma,.) 2 + (Pa,yf) dxdy, (7) 

J 0=1 

where the dimensional constant J7 has a meaning of exchange integral. 

Obviously, equations © and ((SJ) coincide. It is easy to show, that ([5]) coin- 
cides with the equation of motion on a maximal degenerate orbit. Recall, that 
([5]) is obtained when correlations between fluctuations of {P%} are neglected. 
Presumably, equations of motion on other orbits are derived from (|4|) via the 
mean field averaging with more complicate correlation rules. 

A generic orbit of SU(3) is determined by the following equations: 

SabHaVb = m 2 + q 2 
dabcfJ-a^b^c = <?(3m 2 - q 2 ). 

The SU(3)-invariant Hamiltonian on this orbit has the form 

a=l 

- 2\/3q(3m 2 - q 2 )Ha,xr)a^, (8) 

where r\ a is a quadratic form in {/i a }: r\ a = y/§ dabc/J-bl^c- For more details see [5J. 

The Hamiltonian systems on coadjoint orbits of SU(3) serve as classical effec- 
tive models for the spin system of s 1 with biquadratic exchange. Evidently, 
these models describe large-scale (or slow) fluctuations of the mean field. In 
this paper we suppose that the order parameters m, q are fixed numbers. But 
generally speaking, they depend on a temperature T and the interaction con- 
stant J. Taking into account these dependencies, one can consider small-scale 
(or quick) fluctuations of the mean field. 

2.4 Geometrical properties of effective Hamitonians 

Each coadjoint orbit of a semisimple Lie group is a homogeneous space that ad- 
mits a Kahlerian structure. Thus one can introduce a complex parameterization 
of an orbit. For this purpose we use a generalized stereographic projection (for 
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more details see [5]). In the case of group SU(3), the projection is represented 
by the following formulas: 

Ma = - (a + m£ a , Va = ^^ft^S Cq + 2m ^ a , 

W 2 + W 2 + U>3 + WS _ (1 — M>l)(«>3 — «>1«)2) + (1 — Wx)(W3 — WlW 2 ) 

_ . 1"3 — fi>3 — ™2 + W 2 _ . (1 + lBl)(l'3 — IUli02) — (1 + Wl)(w3 — W 1 W 2 ) 

2 ~ 1 V2(l + \w 2 \ 2 + |u> 3 | 2 ) \/2(l + |«Ji| 2 + \w 3 - u>iu> 2 | 2 ) 

. I^2| 2 -|^3| 2 , l-|»l| 2 

w - ; , r- 15 , 1... 19 S3 



1 + |t«2| 2 + |lU3 I 2 1 + l™l| 2 + 1^3 — W 1 W 2 \ 2 

^ _ _ w 2 w 3 - w 2 w 3 5 — I Wl — Wl ^ 

1 + |l"2| 2 + |^«3| 2 1 + |wi| 2 + l™3 — W 1 W 2 \ 2 

_ W 2 + W 2 — W 3 — W 3 _ (1 + M>l)(«>3 — W1W2) + (1 + U>l)(«>3 — W 1 W 2 ) 

~ ^2(l + \w 2 \ 2 + \w 3 \ 2 ) 65 ~ ^(l + \ Wl \ 2 + \w 3 - Wl w 2 \ 2 ) 

_ . w 2 — w 2 + w 3 — w 3 _ . (1 — tDi)(tt> 3 — 10111)2) — (1 — l"l)(lS3 — tDitD 2 ) 

6 ~ ' V2(l + |u) 2 | 2 + |«>3| 2 ) \/2(l + \W!\ 2 + |tw 3 - -u)itu 2 | 2 ) 

_ M>2«>3 + W 2 tU3 ^ _ Wl + «>1 

1 + |lU 2 | 2 + |l"3 | 2 1 + | 2 + \ll>3 — WiW 2 \ 2 

^ 2 - |tu 2 | 2 - |m 3 | 2 ^ 1 + |t»i| 2 - 2\w 3 - mim 2 | 2 

C8 " V3(l + |» 2 | 2 + l«- 3 | 2 ) ~ \/3(l + K| 2 + \w 3 - w lW2 \ 2 ) ' 

The coordinates {w\, W2, W3} (Bruhat coordinates according to [7]) parameter- 
ize a generic orbit of SU(3). In the case of a degenerate orbit, one has to assign 
m = and w\=0, or m = v3 g and W2 = 0. 

In terms of {u> a } the effective Hamiltonians ([7]) and JS| have the form 

«-'/5;«(¥S t ¥?)'* a. ^1,2, 3, ,10) 

where z = ce + is a complex coordinate on the plane obtained from the atomic 
lattice after a large-scale limiting process (see Section HT^j) . The tensor g is non- 
degenerate and positively defined, thus it can serve as a metrics on an orbit. 
Its components {g a p} come from fl7j for a degenerate orbit, and from J8|) for a 
generic one. In terms of the auxiliary vector fields {( a } and {£ a } we have 



w 2 =0 



dcgl _ 1 dCa d(a dog2 _ 1 d£a 9£a 

2 ^ \dw a dwp dw a dwp dw a dwp j ' 

Note, that in terms of {w a } the tensor g does not depend on a particular orbit. 

Being a Kahlerian manifold an orbit of SU(3) possesses a Kahlerian potential. 
For this purpose we use a potential $ of the Kirillov-Kostant-Suoriau form: 

$i = ln(l + |wi| 2 + Iws-w^l 2 ), $ 2 = ln(l + |w 2 | 2 + |^ 3 | 2 ), 
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A topological structure of the orbit is characterized by the second cohomology 
group H 2 of dimension 2. That is why there exist two basis 2-forms, for example 
generated by the potentials $1, $2. Each of them defines a topological charge 



O, - _ / V* — ^ k ( dWa ® W P - dWa dw/3 ) dzAdz k - 1 2 
~" 47r ,/ Z—^dw n dwf3 \ dz dz dz dz ) ' 



1 / ' ;2<I> '' 

a,/3 



On a degenerate orbit only one potential is governing, and only one topolo- 
al chargi 
Evidently, 



gical charge exists. Then the expressions for Q k and 'H^,g 6k differ only in a sign. 



H^ ek > 4nJ\Q k \. (11) 

Hence, on a degenerate orbit a minimum ofH e ff is realized if the equality holds, 
that takes place if{w a } are holomorphic or antiholomorphic functions. Here we 
use an idea of Belavin and Polyakov [2] . 

For a generic orbit we define a topological charge by Q = Qi + Q 2 . In order 
to extend inequality (ITT]) to generic orbits with the topological charge Q, we 
construct an effective Hamiltonian by the following formula: 

8 

Utf ^s^^jCik^+^^^ + ftM.,), (12) 

a=l 

d = m 4 + q 4 - mq(q 2 - to 2 ) + 14m 2 q 2 , 
C 2 = § m 2 + q 2 - 4| mq, 
C 3 = 73 mS + 2q3 ~ f m2( ? + 7E mq 2 . 
In terms of {w a } it is reduced to the form (fTX))) with the metrics 
gon _ 1 ( dC - a 9 Ca , d ^ a d ^ a 
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H&>*kJ\Q\, 



Then we get 

and a minimum ofHPJj? is realized if the equality holds, that takes place if {w a } 
are holomorphic or antiholomorphic functions. 



3 Large-scale topological excitations 

Now we construct a particular class of topological excitations that give min- 
imums to the effective Hamiltonians 7^g gl , and 'Hf™ defined by (fT2|) . We 
describe these excitations by holomorphic functions {w a (z)}. Each set {w%(z), 
w 2 (z), ^3(2:)} represents a mean field configuration of the system in question. 

First, we consider a degenerate orbit of SU(3) with to = 0, q — — where 
a nematic state is realized. In order to satisfy the limiting conditions: /13 — > to, 
^8 9j the functions {w a (z)} have to vanish as z — > 00. Let 

w x (z) = 0, w 2 (z)=^ 2 ^, w 3 (z)=^^, a 2 , z 2 , a 3 , z 3 € C, (13) 
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be a large-scale excitation in the magnet in question. The corresponding mean 
field configuration is obtained by substitution of (|T5|) into |J5J). A behavior of 
the components /13 and /is is represented on the Fig. 1. As z — » 00 the values of 
/^3, /ig tend to to, q. 




_ 2_ 
V3 

Fig. 1. Profiles of the mean field components /X3, in the case of 
configuration fl3|) . 



Suppose Z2 = 23. By shifting of the coordinate z one can easily reduce (I13|) 
to the configuration 

wi(z) = 0, w 2 (z)=<f, w 3 (z)=^f, a 2 ,a 3 eC, 

and calculate its topological charge: 



2i ff (a? 



2 



iirJJc {\z\ 2 +a 2 2 +alY 



dzAdz = l. (14) 



In the case of z 2 7^ z 3 the topological charge equals 2. 

It can be interpreted as follows. Each pole of a mean field configuration rep- 
resents a kind of Bclavin-Plyakov soliton. Each soliton gives a unit topological 
charge. Thus, two distinct solitons have the topological charge 2. No continuous 
deformation take a configuration of topological charge 2 to a configuration of 
topological charge 1. If we allow noncontinuous deformation, then two solitons 
can meet at any point and join into one, at the same time an energy is released. 
From (fTT|) . it follows that the released energy equals 47r J per one pole. 

Note, that the energy of configuration (flU)) does not depend on parameters 
of solitons: 02, £2, 13, 23. It means that the excitation can grow (when |<Z2 1 and 
1 03 1 grow) preserving an energy. Such growth immediately leads to destruction 
of an order in the system. 

One can construct a configuration with more than two solitons: 

n m 

wt(z)=0, w 2 {z) = a 2 /\\_(z - z 2 k), w 3 (z) = a 3 / Y\_(z - z 3k ), (15) 

k=l k=l 

here a 2 , a 3 , {z 2 k}^ =1 , and {^3fc}fcLi are fixed complex numbers. If all values 
{ z 2k}k=i and {z3k}kLi are distinct, a topological charge equals n + m. When a 
pole of the function w 2 (z) coincides with a pole of u>3(z), the topological charge 
decreases by 1. But a coincidence of two poles of the same function (for example 
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W2{z)) does not lead to a decrease of the topological charge. It is easy to see, 
that the minimal energy of configuration ()15D equals AkJ • min(n, m). 

Now we consider a generic orbit, where a ferromagnetic state is realized. 
Suppose m = l and 1 — ^^- ln this case, we describe a mean field in the 
magnet by the effective Hamiltonian "Hf™, defined by (fT2|) . Let 

«>i(*) = i^, «*»(*)=i^ M*)=i^, a k ,z k eC, k = 1,2,3. (16) 

be a large-scale excitation of the mean field. A calculation of topological charges 
gives: 1) Qi =3, Q 2 = 2, if ai ^a 2 ^a 3 or ai=a 2 ^a 3 , 2) Qi=2, Q 2 = 2, if 
ai = a 3 7^ «2, 3) Qi =2, Q 2 — 1, if oi = a 2 = a 3 or a\ ^ a 2 = a 3 . 

4 Conclusion and discussion 

Each generalized Heisenberg Hamiltonian with high powers of the exchange in- 
teraction can be reduced to a bilinear form. By a mean field averaging we 
obtain a classical system from the original quantum one. An averaging of the 
Heisenberg equation gives a Landau-Lifshitz like equation for a mean field. Us- 
ing Lie group apparatus, we construct effective Hamiltonians for the classical 
system with SU(3) symmetry. One of them "H^ s is an SU(3)-analogue of the 
Hamiltonian commonly used in theory of magnetism. In addition, we propose 
another one 'Hz.a , which is biquadratic in the mean field. Further, we construct 
examples of topological excitations that give minimums to the Hamiltonians. 
Such excitations can change their shapes and grow preserving an energy. This 
is a probable scenario for the destruction of an ordered state in a 2D magnet at 
nonzero temperature, that agrees with the Mermin- Wagner theorem. 
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